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Abstract 

We investigate the thermal behavior of quarks and antiquarks interacting via 
a temperature-dependent linear potential. The quarks are constituent quarks 
with dynamically generated masses from the background linear cr-model. 

1 Introduction 

The constituent quark model has been rather successful in meson and baryon spec¬ 
troscopy |]I|. Its main ingredients are massive (m ~ 300 MeV) quark constituents 
interacting with a potential, which is linear at large distances and coulombic at 
short distances. A smaller but important role is played by spin-dependent exchange 
interactions which modify the mass spectra. Recent work on the constituent quark 
model emphasizes the dynamical nature of the constituent quark mass. A chiral 
(T-model P may be used to generate the quark mass thereby allowing the important 
coupling of the constituent quark to its pion cloud. By analysing electron scatter¬ 
ing on the proton an estimate 0 of the pion cloud of the constituent quark has 
been obtained. In ref. |@ its effects on baryon spectroscopy have been worked out. 

The linear cr-model coupled to quarks bears some resemblances to the very popular 
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Nambu-Jona-Lasinio model [^, but it treats the 0+ and 0“ mesons as funda¬ 
mental fields. This is an assumption which is true if the system is analysed with 
a low resolution probe or at moderate temperatures. We would like to argue that 
up to a temperature T k. 250, MeV i.e. beyond the chiral transition temperature 
it is still reasonable to separate {a, vf) degrees of freedom from the rest of other 
mesons. Indications in this direction come from lattice simulations p. Treating 
the thermodynamics of the constituent quarks, we allow a deconhnement transition. 
The thermal gluon fluctuations are modeled to lead to a decreasing string tension 
with temperature. We neglect all feedback from the quark dynamics to the gluon 
dynamics. Therefore the decrease of the string tension is given in our model from 
pure gluon QCD. This way the chiral transition temperature may differ from the 
temperature where the string tension goes to zero. Actual lattice simulations show 
a coincidence of the strongly decreasing value of the quark condensate {qq) and an 
increasing expectation value of the Wilson loop operator. For dynamical quarks the 
Polyakov line operator, however, is no longer an order parameter. The problem of 
double counting the scalar and pseudoscalar degrees of freedom (a, if) as qq bound 
states and as explicit degrees of freedom can be solved by excluding mesons with 
quantum numbers and ^Pq) in the summation over bound qq states. 

The partition function of the constituent quark model is given as: 

Z = J VijVi>V(j)e--fo ( 1 ) 

with the Euclidean Lagrangian 

C = - g{a+ + {^T^)a;V{x 

( 2 ) 

where U{(j)) = ^(00)^. 

Here 'ifj = is the spinor of u and d quarks. The meson helds are combined 
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into the four-component vector 


0 = 


( 3 ) 


The parameters of the linear cr-model at T = 0 are chosen in such a way that the 


minimum of the potential U{(j)) lies at (cxo) = 0.093 GeV. The light constituent 


quarks have a mass of m = 300 MeV and the a mass = 2m. Then the coup¬ 
lings are /Iq = (0.495 GeV)^,A = 2S.33,g = 3.23. These couplings are taken as 
temperature-independent. The potential V between qq in a color singlet state is 
given by Vqq{x — y) = k,\x — y\ — 2a/k with the string tension k = kq = (447.2MeV)^ 
at zero temperature. The second term takes into account the self energy correction in 
the constituent quark model We do not explicitly consider the gluon dynamics, 
but we consider the string tension as a function of temperature. Lattice simulations 
for SUScoior 0 0 have determined the critical temperature Tc = 260 MeV but not 
yet the critical behavior, how the string tension decreases with temperature. For 
large space time dimensions d it is possible to expand around a nontrivial stationary 
point of the string action. In a 1/d -expansion the decrease of the string tension 
has been calculated as a function of temperature |^. The effective string tension at 
hnite temperature is 



( 4 ) 


with 



( 5 ) 


Putting d = 4 and neglecting the higher order corrections would give a transition 
temperature of 220 MeV. We choose the value Tf. = 260 MeV from lattice simu¬ 
lations. The decrease of Ke// at small temperatures is analogous to the universal 
Goulomb correction in addition to the linear potential nr. 

The above simplihed partition function is still rather hard to solve, since the 
bound states formed at low temperatures will grow in size and overlap with increas- 
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ing temperature. The resulting problem of a color correlated system of overlapping 
quark antiquark clusters can probably only be solved with a Monte Carlo calcula¬ 
tion 0. We will limit ourselves in this paper to a simplihed analytical calculation, 
where in the low temperature limit we evaluate the qq partition function as a sum 
over meson bound states. In the high temperature limit we approximate the quarks 
and antiquarks as a correlated gas with the string potential acting as a perturbation. 

The outline of the paper is as follows. In section 2 we calculate the spectrum 
of a single bound qq meson coupled to a heat bath, in section 3 and 4 we compute 
the partition function of the many-quark and antiquark system in the low and high 
temperature schemes given above. In section 5 we address the chiral phase transition 
in mean field theory. 


2 Spectrum of a single qq bound state coupled to 
a heat bath 


In a hrst step we neglect the “elementary” mesons {a, if) and take all mesons as 
bound states of quarks and antiquarks with fixed quark masses. This approximation 
destroys the nice low temperature behavior of the chiral a-model, which is governed 
by the low lying Goldstone pions, but we use this simplification to show our new 
methods to solve the relativistic bound state problem. This technical development 
is one of the main new results of this paper. The same method can also be included 
in more evolved Monte Carlo calculations at finite temperature and finite density. 
Note that the calculation has to be relativistic since with increasing temperature 
the mean field mass of the quarks will decrease, so we cannot use the framework 
of the nonrelativistic constituent model. Instead we propose to treat the partition 
function of relativistic quarks with the help of auxiliary variables. We start with 
the example of a single meson composed of a quark and antiquark with total spin 
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5 = 0 coupled to a heat bath. The quarks have hxed masses m; then 


Zqq = trexp{-/3 + + V^(|rg - rj))} . (6) 


The trick is to rewrite the relativistic Boltzmann factor as follows: 



In the appendix we outline the derivation of the full partition function using the 
same parametrization of the exponential integral as above. Then Zqq is given as 

^ roc poo 

Zqq = tr— / dui / da 2 e~^ 
n Jo Jo 

ii = + /^2 + (^) 

4/ii 4/^2 

This allows to separate relative and c.m. motion in the relativistic two-body prob¬ 
lem. Of course, we still have a static potential without retardation effects in the 
interaction. We also do not know how to include relativistic spin-spin and spin-orbit 
interactions. 

The relative and c.m. coordinates are: 


R = 


r = rq-rq 


fxjfq + fX^fg 

hi + hi 


p = 


PlPq - PlVq 
hi + hi 


P =Pq+ Pq. 


(9) 


With these coordinates the pseudo Hamiltonian Ji. has the form 


TV _ 2,2, 

H — -f /i2 H-7- 


+pV{r). 


2^2 / 1 ^ 

^ H-: 

hi h2 


— + 7:2 ) + y 


p2 


+ 




hi + hi hihi 


(h? + hi) 


( 10 ) 


At this stage it is advantageous to convert to new variables 

h? + hi j 

X = — 7 ^— and y = 


V2f3 


hi - hi 

72/9 


( 11 ) 
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which give 


Zqq 

h 


h 


cm 


h 


rel 


^2 1 
tr— f3 / dx dy ^exp(-/3h) 

TT Jo J-X ~ y 

2 2 

+ 2 ^/\x-y) + 2 ^l\x + y) + 

p2 


25/2a; 

-|- 2 ^eflf • 

TTlred 


( 12 ) 


We see that the separation into a free c.m. Hamiltonian hem and a Hamiltonian of 
relative motion hj-ei has been achieved. The c.m. kinetic energy contains an effective 
“mass” proportional to x, whereas the energy of relative motion is proportional to 
the inverse reduced “mass” mred 


rUred = \f2\x - — 

V X , 


(13) 


In the evaluation of the trace we have to sum over all eigenstates of h^m and hrei- 
The eigenvalues of hem are plane waves 


hem\K) - 25 / 22 , 1 -^)- 


(14) 


The eigenvalues of hrei can be solved numerically: 


hre\^nr,lij') 


^rir/ 


TTli-ed 

^nr,£ 


3/ ^eff 


T 2ydieff I 'h?ir-/(^) ^nr,i^ 


TTlj-ed 


( 15 ) 


For the lowest lying states the coefficients an^/ are as follows: ckqo = 2.34, ckoi = 
3.36, aio = 4.09 and 0:02 = 4.25. 

To calculate the equation of state we need the whole spectrum. Instead of 
numerically calculating the whole spectrum, we take for the higher lying qq states 
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the eigenstates of the harmonic oscillator with the gronnd state adjusted to the exact 
solution. The high lying states become important at high temperatures, where this 
approximation is sufficient to get a qualitatively correct picture 


~ (2/1^ + ^ + 2.34) < 


T?7/red 


2 a/ Keff- 


(16) 


In principal this approximation is not necessary, but it allows us to use analytical 
formulas in the following. With the main quantum number n = 2nr + £ the de¬ 
generacy of harmonic oscillator states is g{n) = {n + 1)(| -|- 1), the energies are 


u{n) = {n + 2.34)^ - 2y/K^ and 


^/o roo f+x p 

Zqq = H - [ dx [ dy 


-fi 


2 ^/^x 


K 


71 Jo J-x y/x-^ — y‘^ 


with 


= \/2x + 


m„ 


+ 


m2 


-f- uj{n). 


2‘^l‘^{x — y) 2^/‘^{x + y) 

The c.m. motion can be integrated out giving the hnal partition function 

25/4pA 


Zqq — 


TT 


( 5 / 2 ) 


./3-0.5 + 1) (^ + l) dxx^-^ 


r-\-x 


dy 


\/x^ — y^ 






We project out the bound state meson masses at T = 0: 


(17) 


(18) 


(19) 


M(nr,i) i—-\ogZ,A ; 


( 20 ) 


In table 1 we give the mass spectrum of relativistic quark antiquark states. The 
bound states are calculated with hxed constituent masses of 300 MeV for the light 
quarks and 450 MeV for the strange quarks. The numerical errors are estimated 
in the following way: The lower values of the masses are obtained from eq. m 
with only the specihc state in the partition function. The upper values are gotten 
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Table 1: Ground state masses M from a relativistic and nonrelativistic calculation 
for light qq and strange quark ss states (in brackets) The experimental p-meson 
spectrum is given in addition. 



ls(0,0) 

lp(0,1) 

2s(l,0) 

2d(0,2) 

relativistic (MeV) 

nonrelativistic (MeV) 

exp (MeV) 

787 ±7 
(1002 ± 7) 
900 
(1064) 

770 

(1020) 

1205 ± 14 
(1391 ± 15) 
1422.0 
(1517) 
1260 
(1285) 

1494 ± 15 
(1665 + 15) 
1793.1 
(1838) 
1450 
(1680) 

1556 ± 15 
(1724 ± 16) 
1875 
(2190) 
1700 


by taking into account hnite temperature corrections using the free meson partition 
function 

Z = (21) 

with Xt = from eq. (|^) for the particular state we are interested 

in. To compute the upper limits in table 1 we apply eq. (plf) with Z = Zqq. The 
minimal temperatures are 1.5 MeV for the ground state and 3.0 MeV for the excited 
states. The obtained spectrum corresponds approximately to the p meson spectrum 
with p(ls), p(2s), p{2d) and the orbital (£ = 1) excitation ai(lp). The agreement 
of the theoretical energies with experiment is rather good. Especially the low lying 
2s-state is improved by the relativistic Hamiltonian. The splitting between the 2s- 
and Is-state of 900 MeV is lowered to the experimentally observed 700 MeV. Note, 
however, that the spin-dependent interactions are missing. For a constituent quark 
mass of 450 MeV we obtain the equivalent 0-states. Since the physical vector meson 
states are well described, we can explore the dependence of the meson masses on 
the choice of quark masses and string tension. Both of these parameters vary with 
temperature. In Eg. 1 we give the dependence of the ground state energy on the 
constituent mass and on the square root of the string tension. In both cases the 
ground state energy drops. In the limit of small quark masses the meson mass is of 
order of a/kq. For vanishing string tension the meson mass converges towards twice 









the quark mass. Here the numerical error of 10 MeV is comparable to the error 
estimates for the realistic meson masses given in table 1. 
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Figure 1: Ground state meson mass as a function of constituent quark mass 

(squares) and as a function of the square root of the string tension (dashed line). 


3 Nonoverlapping (qq)-inesons at low temperature 


Now, we address the problem to calculate the partition function Z of a system of 
7V-quarks and iV-antiquarks. It can be written as 


Z = trexp(—/3H]v) 

N _ N 

t=i i=i 

1 ^ 

+ 2E^(Gi-%)- (22) 
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As before we are limited in our treatment to nonretarded interactions given by the 
static string tension. In addition we have to sum over all possible arrangements of 
flux tubes connecting qq pairs. To our knowledge this is only possible in a fully 
dynamical way with a Monte Carlo model, as has been proposed by the Ontario 
group 1^. In the following we limit ourselves to low temperatures, where the 
number of mesons is small. Therefore we can safely assume that the individual 
mesons keep their identity and do not interact. Consequently we have a iV-“meson” 
partition function with the Tth quark always correlated with the z’th antiquark 

N N 


^ = Yi i dn2,iexp-[H{i)] 


TT 




(23) 


where for the i’th meson 7Y(i) is the same Hamiltonian as in the single meson problem 

32 






4^2,1 




(24) 


After the introduction of c.m. and relative coordinates we transform to the variables 
Xi and Hi as before and get 

\T ,, 

r+^i dvi 




dxj 




exp[—(3h{i)] 


with 


where 


h{i) = V2xi + 


m 


+ 


m 


2^/‘^{xi-yi) 2^/‘^{xi + yi) 


+ hg, + h 


rel ,2 


(25) 


(26) 



Pf 

hs- 

2^PxP 



^rel,2 

Vi 

mred(H 

mred(f) 

= a /2 (xi 


- + KeSTi - 2y/K^ff-, 

yl] 

Xi / ' 


(27) 
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We again use the harmonic oscillator approximation for hj-ei.i given before 


I ki) 2^1‘^Xi ^ ‘2^1‘^Xi ^ 

hreljil^r,^) ( t ■\ 

VrUredU) / 




Oinr,e,i — (2nr + + 2.3381). (29) 

We consider the integration over the auxiliary parameters Xi and yi with the measure 
, I as a summation over intrinsic degrees of freedom of the meson, i.e. 




then each meson is characterized by the intrinsic parameters |r) = \kriUJn^,l,x,y). 
The meson many-body wave function is a symmetrized state dehned by the occu¬ 
pation numbers of each state |r). The free energy of the mesons has the usual Bose 
gas form, but now including the summation over the Xi and yi coordinates 


^ -exp[-/5£(n,/c,a;,|/)]) 


= EE^?( 


ff n k X y 

V- . .V2 r 


dy V 


71 Jo J-x \/x^ — y^ (Zti) 


j ln(l — exp(—/5e)) 


ff ^ 


E ^ exp[-s(3er] (31) 


where 


e{n,k,x,y) = 


9 2 

m m 


/— Hi Hi 

e. = + v'2x + 
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We follow the standard way to evaluate the free energy. We hrst integrate by parts 
and then rewrite the denominator as a power series before finally integrating the 
energy of relative motion numerically. To complete the discussion of the meson gas, 
we also give the expressions for the entropy and energy density 


P 

s 


u 


-F/V- 

~V 


V 




1 dSr 


Sr — 


pdT 

Ts — p. 


E 

S=1 


exp(— 
I372 


(34) 


The resulting pressure of the meson gas is shown in £g. 2 for two cases (a) for 
the ground state mesons only and (b) for the lowest iV = 15 states. One sees that 
because of the large meson masses the pressure is rather small up to T = 150 MeV. 
Around this temperature the large number of degrees of freedom for iV = 15 makes 
a difference. At a temperature T = 250 MeV a Hagedorn transition occurs, when 
the number of relevant states and consequently the pressure increase tremendously. 
One sees that the treatment of the system as a gas of noninteracting mesons be¬ 
comes unrealistic at high temperatures, since the mean distance between mesons 
becomes smaller than the size of the mesons. The Bose distribution yields the mean 
particle number (N) 


im = vJ2T,('2"'‘)3( 

ff ^ 


n)^-—n 7 r f dxx^^'^ [ dy 


TT' 


5/2 


1 00 1 

(35) 


and 


1/3 


Td = 


Att (-/V) 
3 V 


(36) 


The function in fig. 3 is a strongly decreasing function. At T ~ 100 MeV the 

mean distance between mesons is comparable to the root mean square radius ^(r^) 
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Figure 2: Normalized pressure p/T* as a function of temperature T. The lowest 
pressure is obtained for the conhned mesons when only the ground state mesons are 
considered. A diverging pressure (drawn line) results when mesons with = 15 
main qnantum numbers are taken into acconnt. For comparison the free qq gas 
pressnre (open boxes) and the correlated qq pressure (+) are shown. 


of the gronnd state meson 


\/F) = 


2.33 


1/2 


(37) 


Aronnd this temperatnre the approximation of a qnasifree meson gas fails, the 
mesons start to overlap and interact strongly. We also remark that the expression 
for the entropy leads to a divergent term at T = Tc = 260 MeV. This is irrelevant 
because the assumption of independent mesons already fails at mnch lower temper- 
atnres as discussed above. Clearly, for higher temperatures it becomes more and 
more inefficient to treat the system as a meson gas and a calculation in terms of 
qnark variables is more advantageous. This does not mean that the strings and their 
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Figure 3: The meson groundstate rms radius (squares) as a function of temperature 
and the mean distances between mesons (+) and between quarks and antiquarks 
of opposite color Rgc (dashed line). 

attractive forces have disappeared. The strings will arrange themselves according to 
the quark antiquark distributions, instead of extending between fixed pairs of quarks 
and antiquarks. In £g. 3 one can estimate that near T = 140 MeV the mean dis¬ 
tance Rqc between quarks and antiquarks of opposite color is equal to twice the rms 
meson radius. After this energy the correlated qq system starts to be competitive 
in energy with the meson system. This will be treated in the next section. 

4 Correlated quarks and antiquarks at high tem¬ 
peratures 

When the mean distance between quarks and antiquarks becomes smaller than twice 
the rms radius of the meson state, the mesons will no longer keep their identity, i.e. 
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the strings will rearrange themselves in conhgurations with the smallest free energy. 
An exact calculation probably can only be done with a Monte Carlo simulation. 
Here we would like to present an estimate using the Hamiltonian given before. The 
calculation is based on the nearest neighbour saturation of color forces, a concept 
which has been used in the many quark problem before [|^]. In the high temperature 
limit we consider 7^ as a sum of noninteracting quark and antiquark energies 7Y, and 
a perturbation Tix 


n = Ho + Til 


^0 = n + + 1] V^- + 




'Hi = IK't.sIr, - f,| - 


( 38 ) 

( 39 ) 


99 

//' 


Here the summation over r^, distances is understood under the constraint of a 
minimal total length of strings. In thermal perturbation theory we hrst evaluate the 
free energy related to Hq 


Fo = --lnZ,-, = ^ln(l + 


/9 

24 




¥'7 S'- 


( 40 ) 


The mean number of quarks with a color c equals the mean number of antiquarks 
with opposite color c 

d?p 1 


(A'jo) = {N„} = wl 


(27r)^ eHp + 1 ’ 


= 4i/i:(-i) 


s+1 


d^p 


S = 1 


(2vr) 


1 


27r2 ’ 


s + l 


S=1 


l3eps 




We imagine that each quark is localized in a Wigner Seitz cell of radius Rqc- The 
numerical values of Rqc as a function of temperature are shown in £g. 3. The 
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quasilattice of Wigner size cells also of radius Rqc for the antiquarks may be shifted 
relative to the quarks by a distance \l\ < 2Rqc 


I 47r (Ngc) 

\ 3 V 

Then we can calculate a distribution function for a given 
the nearest antiquark with color c at a distance r away. 



( 42 ) 


quark of color c to hnd 


ff(r) = 



J dfq 9{Rqc -\i + rq\) 

J di 0{2Rqc - \i\)S{f- {fq - fq)). 


(43) 


In £g. 4 we give the scaled distribution function g{r/Rqc){r/Rqc)‘^ and compare it 
to e{l - r/Rqc){r/Rqcf. 

These two curves have approximately the maximum at the same place, but the 
correlation function extends to larger distances. With this numerically calculated 
function we calculate the modified free energy in perturbation theory 

F = Fo + 3Nq,Nq,{ni) ^ (44) 

where (Td/) ■ Worr = f g(r/Rqc)r^dr(Keffr — 2ydtgg). The correlation volume 
Worr = dTT ■ R^q, ■ 0.203 is smaller than the naive volume and the mean (r) = 
/ r^drgir/Rqc)/ f r‘^drg{r/Rqc) is unessentially larger than Rgc 


(r) = 1.03i?,e. (45) 

In fig.2 we show the pressure of the correlated quark antiquark gas times 1/T^ (with 
symbols +). It rises rapidly around T = 150 MeV and overshoots the uncorrelated 
pressure at T ^ 175 MeV. Near T = 145 MeV it already exceeds the still very small 
pressure of the meson gas. The meson contribution is damped by roughly 
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Figure 4: The distribution function g{r/Rqc){r/Rgc)"^ (squares) is compared to the 
uncorrelated theta function multiplied by phase space (line). 

compared to the quark excitations which have Boltzmann weights The aver¬ 

age interaction potential changes sign at T = 175 MeV. Of course this value depends 
critically on the subtraction in the qq potential, which gives good agreement 

with the rho-meson in our model. The inclusion of spin-dependent interactions may 
modify the potential. The rather simplihed calculation presented here conveys the 
correct message. It is not the “deconhnement” temperature Tc = 260 MeV for 
pure glue theory which determines the phase transition in the presence of qnarks, 
but the lower temperature T ~ 145 MeV where the quarks and antiquarks liberate 
themselves from the hxed meson conhgurations. At and above this temperature 
the quarks and antiqnarks are still correlated by strings with a string tension of 
= 145 MeV) =400 MeV, bnt these strings change their partners constantly. 
If this pictnre is correct, the heavy quarkonia would still experience a long range 
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potential in addition to the modified Coulomb potential due to screening. Also a 
heavy quark would find additional light antiquarks nearby at no cost in energy. It 
is in this sense that heavy quarkonia dissolve. The Hagedorn transition of the reso¬ 
nance gas occurring at higher temperature is irrelevant, since it happens outside of 
the range of credibility of the model. It would actually prefer a system of meson re¬ 
sonances instead of the quark antiquark gas, but at this temperature the resonances 
are strongly interacting and the independent meson model has lost all its theoretical 
foundation. 


5 Chiral dynamics of constituent quarks 


Having set up the formalism for a fixed constituent quark mass we are ready to 
consider the variation of the constituent quark mass due to the quark condensate. 
We can use the linear a-model to generate the quark mass m = g{a) where (a) is 
determined from the minimum of the free energies either in the resonance gas phase 
or in the correlated quark antiquark phase plus the potential energy of the mean 
(T-field 




ff ” 


dxx^^'^ f dy 


-X y/x^ — ^ 


E 


^-SpEr 


+U({a}) 


( 46 ) 


where 


Sr — (2,nr T 2.34) 


^lAT) 


+ V2x + 


iaW)) 


\ V2(x - y‘^/x) 


2i/k, 


eff(T} 


+ 


2^/‘^{x — y) 2^/‘^{x + y) 


(47) 


Here one sees how the relativistic bound state calculation takes into account the 
mean field masses g{(7) of the quarks via the integration over the auxiliary parame- 
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ters X and y. 


In the high temperature phase of correlated quarks and antiquarks with variable 
masses the free energy reads 



-^v j lii(l + exp(-l3^p‘ + (g{a}y )) 
+3N,M<y))Nd9W){H,)'^ 


(48) 


where also the interaction energy and the correlation volume depend on the con¬ 
stituent quark mass m = as described in section IV. The mean held quark 

mass is determined by the condition that the free energies F{{a)) are minimal, i.e. 


dF 

W) 


(cr)=(To 


0 ; 


m = gao- 


(49) 


In £g.5 we give the dependence of the constituent quark mass on the temperature for 
quark matter and resonance matter. In both cases we have a first order phase transi¬ 
tion with Tchirai(C'Q) = 147 MeV, Tchirai(7?G) = 196 MeV. Until 120 MeV very little 
change occurs. This is different from calculations of the NJL type with nonconfining 
interactions. In the line of argumentation given before we do not attribute great 
significance to the higher chiral phase transition of resonance meson matter, proba¬ 
bly a more accurate calculation including the low lying tt- and cr-degrees of freedom 
explicitly would make the two chiral transitions coincide. What is very interest¬ 
ing is the almost coincidence between the deconfinement temperature T^ec = 143 
MeV, where the resonating meson gas goes over into correlated quark matter and 
Tchirai{CQ) = 147 MeV. This coincidence is surprising given the simplicity of the 
theoretical calculation. It is dependent on the choice of original constituent quark 
mass and the behavior of the string tension with temperature. It may indicate that 
the constituent quark mass is intimately related to the gluon dynamics and the 
choice of critical behavior of the string tension is perhaps not so far off from reality. 
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Figure 5: The constituent quark mass as a function of the temperature is shown 
for the correlated quark antiquark (lines) and for the resonance gas (squares). 

6 Conclusions 

In the framework of the chiral constituent model we have presented a relativistic 
calculation of the pressure and the effective quark mass of quarks and antiquarks at 
hnite temperature. At low temperature the constituents are bound in mesons which 
change their composition only very slowly. Due to the decreasing string tension the 
quark antiquark bound states have a slightly decreasing mass. In future work we 
would like to investigate further the p meson spectrum as a function of temperature. 
At temperatures around T = 120 MeV the approximation of noninteracting mesons 
becomes inaccurate and a description in terms of liberated quarks and antiquarks 
connected by varying string conhgurations becomes more appropriate. At a temper¬ 
ature T ^ 145 MeV both a deconhnement and chiral transition occurs to massless 
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fermions which prefer arranging their string to the nearest antiquark of opposite 
color nearby. This temperature is distinctly lower than the quark mass transition in 
the resonance gas occuring at 196 MeV. Important fluctuations of the chiral helds 
are still missing in the calculation so far. Especially the Goldstone vf-fields con¬ 
tribute an important part of the pressure at low temperatures. The addition of 
these extra degrees of freedom will be done in another paper. In the resonance gas 
the summation over gg-bound states should exclude the ^S'o(7f) and ^Poicr) states. 
Therefore, for the ground state ^Si and ^S'o state the new degeneracy factor becomes 
g'(0) = 16 — 3 = 13 and g'{l) = 48 — 1 = 47 for the first excited states, where the 
48 states come from the degenerate 4 isospin Pi,^ P 2 and ^Pi channels in the 
harmonic oscillator model. To evaluate the fermion determinant in the presence of 
a and vf fields it is advantageous to use the heat kernel method which has been 
developed for finite temperatures in ref. |^. The methods described in this ar¬ 
ticle are not limited to finite temperature. They can as well be applied to hnite 
baryon density, especially to the nuclear problem which has the most promising fu¬ 
ture. Until now it has resisted any explanation in terms of fundamental quark and 
string degrees of freedom. Constituent quarks, however, have the correct dynamics 
in terms of mesonic chiral interactions and confining string dynamics to make the 
problem accessible. 


Appendix 

The density matrix for the relativistic two body scalar problem with interaction V 
neglecting retardation and backward motion has the form: 


Pqq = + m2 + + m2 + V{\rg - fg]))]. (A.l) 

It has a path integral representation 
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Pj,- = /l’VoI>V^exp|-^ [J\/00 + ^0 +$\/oo + ^j + V(|V^ - Vn|)] fr], 

(A.2) 

This integral can be divided up into a product of integrals over small At inter¬ 
vals which are reformulated with auxiliary helds. The integrals with square roots 
are transformed into integrals without square roots with the help of the following 
integral: 


exp 


a = 


I dxexp-{x^ + ^) 


(A.3) 


Eq. (|A.3|) can be proven by evaluating the integral I and its derivative with respect 
to a 

I = J dxexp[-{x - (A.4) 

(A.5) 


dJ _ 0 

da 


Since the derivative of I with respect to a is zero, the value of I is the norm of the 
Gaussian integral I {a = 0). After a rescaling = fifAr and a suitable choice of 
the measure we obtain the form 


Multiplying the individual time slices with each other one generates the full density 
matrix 


(^u +oo))-r] (A.6) 


exp[-m)Jl + r/Ar] = Af J \ia^exp[-{n^ + 


Pqg = Af' J (A.7) 

exp[-y^ [pi + PI + + Tq ) + ^(l+^<? ) + V{\fq - fq\)]dT. 

In general the auxiliary helds are functions of r. We propose as an approximation to 
keep only auxiliary helds in the functional integral which are independent of time. 
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Then the path integral over can be executed by solving the corresponding 

two-body Schroedinger equation 

^ ^ + '^(1^9 - 1)]^ 5 ) = ’ /^ 2 )'h(fg, fg). (A. 8 ) 

4/ii 4/i2 

The hnal integration over the auxiliary helds becomes a two-dimensional integration 
pqq= W j [/^^ + /^i + ^^ + ^^ + ^(hoo,he)]M. (A.9) 

Finally the partition function can be written in the form of eq. 8 of section 1 by 
the substitution /5/if —/if and /5/if —/if. The normalization factor in front of 
the integral in eq. 8 is adjusted to reproduce the free partition function without 
potential. 
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